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Abstract 

We consider the deformations of Whitham systems including the "dispersion 
terms" and having the form of Dubrovin-Zhang deformations of Frobenius mani- 
folds. The procedure is connected with B.A. Dubrovin problem of deformations of 
Frobenius manifolds corresponding to the Whitham systems of integrable hierar- 
chies. Under some non-degeneracy requirements we suggest a general scheme of the 
deformation of the hyperbolic Whitham systems using the initial non-linear system. 
The general form of the deformed Whitham system coincides with the form of the 
"low-dispersion" asymptotic expansions used by B.A. Dubrovin and Y. Zhang in 
the theory of deformations of Frobenius manifolds. 

1 Introduction. 

The classical Whitham method (Pfl El El E]) is connected with the slow modulations of 
the exact periodic or quasiperiodic solutions of non-linear PDE's 

^ l (V> ¥t, <P X , <Ptt, <Pxt, <Pxx, ■ • •) = , i = l,...,n (1.1) 

where (f = (ip 1 , . . . , (p n ). 

It is assumed that the system (jl.l)) admits the finite-parametric family of exact solu- 
tions 

(p'(x,t) = ¥(k(V)x + w(U)t + O , U) (1.2) 

where 6 = (8 1 , . . . ,6 m ), $ l (0, U) are smooth functions 27r-periodic w.r.t. each 6 a , 
k(U) = (^(U),...,^" 1 ^)), w(U) = (^(U),...,^™^)) are "wave numbers" and 
"frequencies" of the solution, U = (U 1 , . . . , U N ) are parameters of the solution, and 
0q = (#o, . . . , 9™) are arbitrary initial phases. 

The functions <&(0, U) satisfy the nonlinear system 

F l ($, u a (\J)® e «, k p (U)^ 00 , ^(U)a;*(U)* fl ^, . . .) = (1.3) 
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We can introduce the families Ak, w and the full family A = (J Ak,^ of the functions 
$(0, U) satisfying the system ()1.3|) in the space of 27r-periodic w.r.t. each 9 a functions. 
Let us choose (in a smooth way) at every (U 1 , . . . , U ) some function <fr(0, U) as having 
"zero initial phase shifts" and represent the full family of m-phase solutions of system 
(TTTJ) in the form (JO}. 

In Whitham method we make a rescaling X = ex, T = et (e — > 0) of both variables 
x and t and try to find a function 

S(X,T) = (S\X,T),...,S m (X,T)) (1.4) 

and 27r-periodic functions 

¥(0,X,T,e) = ^^ fc) (0,X,T)e fc (1.5) 

k>0 

such that the functions 

<P\0,X,T,e) = &nW±Q- + 0,X,T,e) (1.6) 

satisfy the system 

e0 T , e<l>x, e 2 cj> TT , . . .) = (1.7) 

at every X, T and 0. 

It is easy to see that the function \&( O )(0, X, T) satisfies the system (jl.3|) at every X 
and T with 

fc a = S£ , u a = S% 
and so belongs at every (X, T) to the family A. We can write then 

¥ {o) (0,X,T) = ¥(0 + o (X,T),U(X,T)) 

and introduce the functions U U (X, T), 9q(X, T) as the parameters characterizing the main 
term in (jl.5[) which should satisfy the condition 

[fc«(U)] T = [uf{U)] x (1.8) 
The functions X, T) are defined from the linear system 

L^ j {1) (0,X,T) = f (l) (0,X,T) (1.9) 

where 



dF i 

L ) = L U, T)j = w (* (0) (e,x,T),...) + 
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+ — (* (o) (0, X, T), . . .) W «(X, T) — + ^ (* (0) (0, X, T), . . .) ^(X, T) — + . . . 

(1.10) 

is the linearization of system (jl.3j) and f(i)(0, X, T) is discrepancy given by the expression 



f {1) (0,X,T) = - (vfr (o) (0,X,T), ...) ^ o)r (0,X,T)- 



- — (* (0) (0, X, T), . . .) 9> m (0, X T) - 



* (o) (0,X,T), ...) (2^(X,T)M/J 0)eQT + ^(X,T)^ 0)e/3 ) - ... (1.11) 



We have here 

d d d d d d 

of ~ ~du^ + u ^ T m^ ' 9x " Ux mF + U ^ X W 

for the functions 

*\ 0) (0,X,T) = $ l (0 + o (X,T),U(X,T)) 

We will assume that fc a and co> a can be considered (locally) as the independent param- 
eters on the family A and the total family of solutions of (jl.3|) depends (for generic k a , 
u a ) on N = 2m + s, (s > 0) parameters U u and m initial phases 9q. 

Easy to see that the functions &g a (6 + 6 (X,T),XJ(X,T)) and 
V^^ga(0 + 6o(X, T), U(X, T)) where £ is any vector in the space of parameters U u 
tangential to the surface k = const, uj = const belong to the kernel of the operator 

Let us put now some "regularity" conditions on the family (J1.2)) of quasiperiodic 
solutions of (jl.lj) 

Definition 1.1. 

We call the family J 1.2]) the full regular family of m-phase solutions of M.l)) if: 

1) The functions <Ev(#,U) , $[/^(0,U) are linearly independent (almost everywhere) 
on the set A; 

2) The m + s linearly independent functions $£)«(#, U) ; V^$(0,U) (V^k = ; 
V^cj = 0) give the full kernel of the operator LLj, . (here 6q = 0) for generic k and uj. 

3) There are exactly m + s linearly independent "left eigen-vectors" q = 

l,...,m + s of the operator L 1 ^. (for generic k and u) corresponding to zero eigen 
values i.e. 

»2tt r 2-K 



pZtt pZtt r] Ta 
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for any periodic ipi(0). 

It's not difficult to see that the Definition 1.1 is connected with the regularity properties 
of the sub-manifold A given by the set of functions <&(0, U) in the space of 27r-periodic 
functions. In fact, for our purposes we can use also the weaker definition of the full regular 
family of m-phase solutions of (Jl.lj) . Namely, let us represent the space of parameters 
U in the form U = (k, u>,n) where k are the wave numbers, uj are the frequencies of 
m-phase solutions, and n = (n 1 , . . . ,n s ) are some additional parameters (if they exist). 
Let us give now the "weak" definition of the full regular family of m-phase solutions in 
the form: 

Definition 1.1'. 

We call the family A the full regular family of m-phase solutions of M-l\) if 

1) The functions k, uj, n), & n i(6, k, uj, n) are linearly independent and give (for 
generic k and uj ) the full basis in the kernel of the operator L 1 .^ k w 

2) The operator L^ 6o k w n j has (for generic k and uj ) exactly m+s linearly independent 
"left eigen vectors" 

depending on the parameters U in a smooth way and corresponding to zero eigen-values. 

We will assume now that the system (jl.lj) has a full regular family of quasiperiodic 
m-phase solutions in strong or weak sense in all our considerations. 

To find the function ^/^(6,X,T) we have to put now the m + s conditions of or- 
thogonality of the discrepancy i^(0,X,T) to the functions k [u(xt)](^ + ^o(X, T)) 

/2tt i-2-k rj m ft 
■■■jf ^CT.( + «( I ' T ))/(1)( ' I ' T )(^ = ( L12 ) 

The system (11.12)1 together with p. 8)1 gives m + (m + s) = 2m + s = N conditions 
at each X and T on the parameters of zero approximation Vl/(o)(0, X, T) necessary for the 
construction of the first e-term in the solution ()1.5)) . 

Let us prove now the following Lemma about the orthogonality conditions 1)1.12)1 : 

Lemma 1.1. 

Under all the assumptions of regularity formulated above the orthogonality conditions 
\1.12}) do not contain the functions 9q(X, T) and give just the restrictions on the functions 
U V {X, T) having the form 

Cl q) {U)U£ - D^\V)U V X = (1.13) 
(with some functions Cj> q) (U), D { u q \\J)). 
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Proof. 

Let us write down the part of the function f(i) which contains the derivatives 
eSr(X,T) and 6^ X (X,T). We have from (fl~TTl) 

X, T) = - — (* (0) , . . .) < - (* (0) , . . .) C - 

~ ITT (*(°)> • • •) 2 ^ T ) ^ ~ FT (*(P)> • • •) 2 ^ T ) *(o)^ & - • • • 

Let us choose again the set of parameters U in the form 

U = {k\ . . . ,k m ,u\ . . . ,uj m ,n\ . . . ,n s ) 

where k a are the wave numbers, uj a are the frequencies of m-phase solutions and (n 1 , . . . , n s ) 
are additional parameters (except the initial phases). 
We can write then 



V + 



+ 



Vox 



The derivatives dF % /duj 13 and OF 1 /dk 13 are identically zero on A according to (jl.Hj) . 
We have then 



2tt 



2- 



d m e 



= o 



since all K^ q \&, X, T) are the left eigen- vectors of L with zero eigen- values. 

It is easy to see also that all 0q(X,T) in the arguments of $ and will disappear 
after the integration w.r.t. so we get the statement of the Lemma. 

Lemma 1.1 is proved. 



The system 



dk a 



X ) 



a = 1 m 



C®{\S)T% = D^(V)U X , q=l,...,m + s (1.14) 

is called the Whitham system for the m-phase solutions of system (jl.lj) . 

Let us note that we have rank\\dk a / dU v \ \ = m according to our assumption above. 
In the generic case the derivatives U!f can be expressed through U x and the Whitham 
system (II . 14)1 can be written in the form 
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U$ = V»(V)U X ,1/,//= 1,...,N (1.15) 

where V^([7) is some N x N matrix depending on the variables U l , . . . , U N . 

Let us say that quite often the system (jl.l|) can be written in the evolution form 

ft = QX<P, <Px, Vxcc, •••) (1.16) 

For systems ()1.16|) the form ()1.15|) of the corresponding Whitham system has then a 
natural motivation. 

We will assume here that if the conditions ()1.14j) are satisfied then the system (jl.9j) is 
resolvable on the space of 27r-periodic w.r.t. each 6 a functions. The solution Vl/(i)(0, X, T) 
is defined then modulo a linear combination of the "right eigen-functions" of L\ XT ^ 
(\P(o)0«, $( )„i) introduced above. According to common approach ([HIE]) we can try 
to use the corresponding coefficients to make the systems analogous to (jl.9j) resolvable in 
the next orders and try to find recursively all the terms of series ()1.5|) . 

Different aspects and numerous applications of the Whitham method were studied in 
many different works (pQ-jHI]) 1 and the Whitham method is considered now as one of the 
classical methods of investigation of non-linear systems. 

It was pointed out by G. Whitham ([UEIIS]) that the Whitham system ()1.14j) has a 
local Lagrangian structure in case when the initial system has a local Lagrangian structure 



5 J J C((p, (fi t , cp x , . . .) dxdt = 



on the space {(fi(x, t)}. 

The procedure of construction of Lagrangian formalism for the Whitham system (jl,14j) 
is given by the averaging of the Lagrangian function C on the family of m-phase solutions 
of system (jl.lj) ([H 12 IS])- Let us note also that in the case of presence of additional 
parameters n l the additional method of Whitham pseudo-phases should be used. 

The important procedure of averaging of local field-theoretical Hamiltonian structures 
was suggested by B.A. Dubrovin and S.P. Novikov ([U UHl EH EI])- The Dubrovin- 
Novikov procedure gives the local field-theoretical Hamiltonian formalism for the Whitham 
system (|1.15J) in the case when the initial system ()1.16j) has a local Hamiltonian formalism 
of general type. The Dubrovin- Novikov bracket for the Whitham system has a general 
form 

{U"(X), U»{Y)} = g>" M (U) 8'(X — Y) + 6^(U) U x 8(X — Y) (1.17) 

and was called the local Poisson bracket of Hydrodynamic type. The theory of the brackets 
(|1.17j) is closely related with differential geometry ([TJJ 123 Ell) an< l i s connected with 
different coordinate systems in the (pseudo) Euclidean spaces. Let us say also that during 

1 We apologize for the impossibility to give here the full list of numerous works on Whitham method. 
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the last years the important weakly-nonlocal generalizations of Dubrovin-Novikov brackets 
(Mokhov-Ferapontov bracket and Ferapontov brackets) were introduced and studied 

E3E3E3EH1EIIEH1). 

The Hamiltonian structure (|1.17|) for the systems (|1.15j) has a direct relation to the 
integrability of the systems of this class. Thus it was conjectured by S.P. Novikov that any 
diagonalizable system (|1.15|) Hamiltonian with respect to the bracket of Hydrodynamic 
Type is integrable. The conjecture of S.P. Novikov was proved by S.P. Tsarev (|51j) 
who suggested the "generalized Hodograph method" for solving the diagonal Hamiltonian 
systems (I1.15|) . Let us say that the Tsarev method has become especially important for 
the Whitham systems corresponding to integrable hierarchies and provided a lot of very 
important solutions for such systems in different cases. We note here, that the Whitham 
systems of integrable hierarchies can usually be written in diagonal form ( El 1213) and 
admit the (multi-) Hamiltonian structures given by the averaging of the Lagrangian or 
the field-theoretical Hamiltonian structures of the initial system. 

During the last years the theory of compatible Poisson brackets fll.l7j) and their defor- 
mations in connection with Quantum Field Theory was intensively developed ( [oDl I?H1 
E21E3ISIESIE9E3ESIE21IZDIIZI1)- Namely, the theory of compatible Poisson brackets 
(I1.17|) plays the main role in the theory of "Frobenius manifolds" constructed by B.A. 
Dubrovin and connected with the classification of the Topological Quantum Field Theo- 
ries (based on the Witten - Dijkgraaf - Verlinde - Verlinde equation). Every "Frobenius 
manifold" is connected also with the integrable hierarchy of Hydrodynamic Type 

K = Vg )M (U)l/£ (1.18) 

having the bihamiltonian structure in Dubrovin-Novikov sense with a pair of Poisson 
brackets ([TT7]) . 2 

The " e-deformations " of the integrable hierarchies and the Poisson brackets (jl,17|) 
are connected with the "higher genus" corrections in the Quantum Field Theories and 
are being intensively investigated now ([HB1 EH EE EH ETJ EE EH UHl El)- According 
to B.A. Dubrovin and Y. Zhang the e-deformation of Frobenius manifold is given by 
the infinite series polynomial w.r.t. derivatives Uj, Uxx, Uxxx, • • • and represent- 
ing the "small dispersion" deformation of the hierarchy (jl,18|) and the corresponding 
bihamiltonian structure of Hydrodynamic Type. The deformed hierarchy and bihamilto- 
nian structure should possess also the additional properties of Frobenius manifolds ("unit 
vector field", Euler field, r-symmetry etc. ...) and give then the "deformation" of the 
corresponding Frobenius manifold Q631 EI3 ETJ)- The general problem of classification of 
deformations of bihamiltonian hierarchies (jl,18|) is being also intensively studied by now 
and very important results were obtained recently in this area IoTH 1701 EH)- 

The dispersive corrections to the Whitham systems were first considered by M.Y. 
Ablowitz and D.J. Benney ([5], also |E]-[Z1) where the first consideration of multi-phase 

2 The hierarchy Ijl.f 8|) and the corresponding pair of Poisson brackets possess also some additional 
properties (the existence of "unit vector field", Euler field, T-symmetry etc. ...). 
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Whitham method was also made. As was shown in j3] the dispersive corrections to 
Whitham systems can naturally arise and, besides that, can be generalized to the multi- 
phase situation. 

In this paper we will consider the deformations of the Whitham systems (|1.15|) having 
Dubrovin-Zhang form, i.e. the dispersive corrections to (|1.15|) containing the higher X- 
derivatives of the parameters U. The problem in this form is connected with the problem 
set by B.A. Dubrovin which is formulated as the problem of deformations of Frobenius 
manifolds corresponding to the Whitham systems of integrable hierarchies. B.A. Dubrovin 
problem contains both the problem of deformation of Whitham systems (jl.15)) and the 
corresponding bi-Hamiltonian structures of Hydrodynamic Type giving the dispersive 
corrections to these structures. 

We will consider here the first part of B.A. Dubrovin problem and suggest a general 
scheme of recursive construction of terms of the deformation of the Whitham system (jl,15j) 
using the initial system or (jl,16|) . The second part of B.A. Dubrovin problem will 
then be considered in the next paper. We will not require, however, the " integrability" of 
the system (11.15)) and all the considerations below can be applicable for both integrable 
and non-integrable systems (jl.l|) . f)1.16j) . Thus, we do not require also the bihamiltonian 
property of the Whitham system ()1.15j) . 

In the next chapter we will make more detailed investigation of the asymptotic series 
(I1.5J) and describe the construction of the deformation procedure for general Whitham 
systems. 

2 The deformation of the Whitham systems. 

Let us start again with the description of the construction of asymptotic series (|1.5j) 
connected with the Whitham method. We will assume that we have the initial system 
having the general form ()1.1|) which has an A-parametric (modulo the initial phases 9$ ) 
family of m-phase solutions 

<p\x,t) = &{k(U)x + w(U)t + O , U) (2.1) 
with the functions <&(6, U) satisfying the system 

Lu a (XJ)<S> e «, k p (U)^ p, ^(V)co 8 (V)^ e s, ...) = (2.2) 

(and 27r-periodic w.r.t. all 9 a ). 

The choice of the functions $*(#, U) at each U is defined modulo the initial phase 6 
and the full family of solutions of (|2.2|) is given by all the functions $(0 + 6 , U) with 
arbitrary Q = (9^, . . . , 9™). We will just assume here that the choice of $*(0, U) is 
smooth on the family of m-phase solutions A. 

We will assume also that the parameters (k, uj) are independent on the family A such 
that N > 2m. It will be convenient to use the parameters (A; 1 , . . . , k m , u , . . . , u m , n 1 , . . . , n s ) 
on the family A where k a are the "wave numbers", u a are frequencies, and (n 1 , . . . ,n s ) 
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are additional parameters (if they present). The linearly independent solutions of the 
linearized system ()2.2j) will then be given by m + s functions &g a (0, U), 3> n i(0, U). 

As we already said above we assume that the family A represents the "full regular 
family" of m-phase solutions of (|2.3|) such that all the requirements (l)-(3) formulated 
in the Definition 1.1 are satisfied. We have thus exactly m + s functions <&g a (6, U), 
<& n i(0, U) giving the basis in the kernel of linearized system (|2.2j) and exactly m + s 
linearly independent functions k^JO + 6 ) giving the basis in the kernel of the adjoint 
operator (for generic k, uj) and depending in a smooth way on the parameters (k, u;,n). 

As we mentioned already we make a rescaling of coordinates X = ex, T = et and 
try to find the solutions of system 

F\cp, e<p T , eipx, e 2 <p TT ,e 2 <p XT ,e 2 <p X x, ■ ■ ■) = , i = l,...,n (2.3) 
having the form 

<F{B,X,T,e) = &(?&>?1 + 0,X,T,e) (2.4) 



¥(0,X,T,e) = X> l (fc) (0,X,T)e fe (2.5) 



k>0 



The function ^^(6, X, T) belongs to the family A at every fixed X and T and the 
compatibility conditions (|1.12j) for the system 

L)* j {1) (0,X,T) = r {l) (0,X,T) (2.6) 

give the Whitham system (|1.13|) on the parameters U(X, T) of the zero approximation. 
The function ^f^(6, X, T) is defined from the system (J2.6)) modulo the linear combination 
of functions $>gn(0, U), <& n i(6, U) at every X and T. 

All the other approximations \& X, T) satisfy the linear systems 

L)¥ ik) (e,X,T) = f( k) (0,X,T) (2.7) 
where the functions f(^) (0, X, T) represent the higher order discrepancies given by system 

(Q. 

The compatibility conditions of the systems (|2.7j) (k > 2) give the restrictions on the 
"initial phases" 8q of the zero approximation and on the previous corrections ^^'){6, X, T). 

Let us make now one more general assumption about the systems ()2.6|) and ()2.7|) . 
Namely we omit here the special investigation of the solvability of systems ()2.6|) . ()2.7|) on 
the space of 27r-periodic functions and assume that the orthogonality conditions 

-jf «mx,T)]i(0 + <>o(X,T)) fl k) {0,X,T) — = (2.8) 

give the necessary and sufficient conditions of solvability of these systems. So, we will 
assume that under the conditions ()2.8|) we can always find a smooth 27r-periodic solution 
of (|2.6j) . (J2.7|) defined modulo the linear combination 
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5>f fc) (X,T) & e «(0,X,T) + J2 d {k)(X,T) $ n i(0, X, T) (2.9) 

a=l 1=1 

Let us consider now the construction of the asymptotic series (12. 5 j) . We note first that 
the solution ^^(0, X, T) is defined from the system ()2.7j) modulo the linear combination 
(12.91) which is equivalent in the main order to the addition of the values e k+1 c?U (X, T) 
to the phases S a (X, T) of the zero approximation ^^(0, X, T) and to the addition of 
the values e k d 1 ^ (X, T) to the parameters n l (X,T) of ^^(0, X,T) according to the for- 
mulae ()2.4j) . (|2.5J) . Easy to see also that we can add the values e h+1 c" k ^ x (X,T) and 
e k+1 c" fc ^ T (X, T) to the parameters k a and u) a in the U-dependence of \I/(o)(0, X, T) which 
does not affect the k-th order of e in the series (12.51). 



Let us change now the procedure of construction of series (J2.5)) in the following way: 

1) At every step k we choose the solution ^^(0, X, T) in arbitrary way; 

2) We allow the regular e-dependence 

S a (X,T,e) = J2 S U X > T ) ek 



k>0 

It 



l (X,T,e) = 5>; fc) (X,T)e A 



n 

k>0 



of the phases S a and the parameters (k, uj,n) of the zero approximation ^/^(6,X, T) 
such that 

k a (X,T,e) = S$(X,T,e) , u a (X,T,e) = S%(X,T,e) 

3) We use the higher orders Sfy, n\ k ^ of the functions S a (X, T, e), n l (X, T, e) to provide 
the orthogonality conditions ()2.8|) for the systems ()2.7|) . 



We can write now the asymptotic solution of ()2.3j) in the form 



[0,X,T,e) = $M S(X,T,6) + 0, S x (X,T,e), S T (X,T,e), n(X, T, e) ) + 

S(X,T,e 



k>l 



+ £*U^^ + ^ T efc ( 2 - 10 ) 



Let us note that the phase shift 0o{X,T) of the initial approximation becomes now 
the e-term in e-expansion of the function S(X,T,e). 



10 



It's not difficult to see that the series (|2.10|) can be always represented in the form 
()2.4|) - ([2.5|) and give the same set of the asymptotic solutions of system ()2.3jl . However, the 
series (J2.10|) contains a big " renormalization freedom" since we can change in arbitrary 
way the higher orders corrections S( 2 )(X,T), S( 3 )(X,T), . . ., n^(X,T), n( 2 )(X,T), ... 
and then adjust the functions *&(k) in the appropriate way. 

Let us now "fix the normalization" of the solution (J2.1(Jj) in the following way: 



We require that the first term 
'S(X,T,e 



+ 0, S x (X,T,e), S T (X,T,e), n(X,T,e] 



of (|2.1()jl gives "the best approximation" of the asymptotic solution (|2.1()|) by the modu- 
lated m-phase solution of (jl.lj) . To be more precise, we require that the rest of the series 
(glint is orthogonal (at every X and T) to the vectors $^(0, X, T), $ n ;(0, X, T) 
"tangent" to A at the "point" [S(X, T, e), n(X, T, e)]. 

Thus we will put now the conditions 



... ^(0, 8 X , St, n) 

Jo Jo i=l 



J2 %)(0, X, T) e< 



k=i 



d m 9 
(2vr) m 



(2.11) 



p27T p27T 71 

... ^$i,(0,S X) S T) n) 
Jo Jo i=1 



J2 %)(0, X, T) e* 



fc=i 



d m e 

(2vr) m 



(2.12) 



(for all e). 



As we saw above, the functions vE^m (0, X, T) are defined modulo the linear com- 
binations of the functions 3v(0, S( )x? &(o)Ti n (o)) ; 3v(^ &{0)x, S(o)t, n (o)) which are 
linearly independent for the full regular family of m-phase solutions of (jl.lj) . Using this 
fact it's not difficult to prove that the conditions (j2.11j) - (j2.12j) fix uniquely all the terms 
S(fc), nm, and 4*^(0, X, T) (k > 0) for a given solution (|2.1(J|) . 

Let us say now that the choice of the normalization (j2.11j) - (j2.12j) is not unique. In 
particular, it depends on the choice of the variables ip l (x,t) for the "vector" (i > 1) 
systems (II. lj) . We will speak about the normalization (|2.11|) - ()2.12|) as about one possible 
way to fix the functions 4f^(0, X, T). 



The solution <fi (0, X, T) can now be defined as the asymptotic solution of the system 
(j2.3j) having the form ()2.10j) which satisfies the conditions 



r-2n 



r2ir n / 

\ J2% a (0, S x , S T ,n) pi 
Jo i=1 V 



- -, X, T, e 



d m e 



(2 
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Jo Jo i=1 



0-|, X,T,e) - &(0,S x ,S T ,n) 



d m 9 



(for all e). 



To find a solution in the form (j2.10|) we substitute the series (j2.10|) in the system (|2.3j) 
and try to find the functions ^i k J0, X, T) from the linear systems 

%o),n ( o)]^)(^ r ) = fl)(^X,T) (2.13) 

analogous to (|2.7|) . 

The functions fj k ^(0, X,T) are different from the functions f^(0, X,T) since we in- 
cluded "a part of e-dependence" in the first term of ()2.10|) . The functions i^)(0,X, T) 
should be orthogonal to the functions ,(0,XT), q = 1, . . . , m + s and we use the 

° l°(0)i n (0)J v ' ' ' * 

conditions (|2.11j) - (J2.12)) for the recurrent determination of \& (&)(#, X, T). The functions 
S^lX, T), n^(X, T) are used now to provide the resolvability of the systems (|2.13J) in 
all orders of e. 

Let us investigate now the systems arising on the functions S( )(X, T), S(i)(X, T), 
. . ., ii(o) (X, T), ri(i) (X, T), .... Let us note that we have in our notations S(X, T) = 
S(o)(X, T), 0o (X, T) = S(i)(X, T). We saw in Lemma 1.1 that the function 0q(X, T) does 
not appear in the solvability conditions of the system (jl.9|) . For the asymptotic solution 
written in the form ()2.10j) with the normalization conditions fl2.11|) - ()2.12|) we can prove 
here even stronger statement. 

Lemma 2.1 

The functions S^ k JX,T), n, k JX,T) do not appear in the expression for the discrep- 
ancy f(fc)(0, X, T) and do not affect the solution *S?(ty(0,X,T). 

Proof. 

The way how we get the discrepancy i^(0, X,T) can be described as follows: 

1) We substitute the solution (|2.1(Jj) in the system (|2.3|) : 

2) After "making all differentiations" we can omit the argument shift S(X, T, e)/e in 
all functions depending on 0; 

3) Then collecting together all the terms of order e k we get the system ()2.13|) . 

It's not difficult to see that in this approach the k-th order of e containing the functions 
S(fc)(X, T) or n(fc)(X, T) have the form 

dF l Qa dF* dF' l 

dk^ b ^ x ' d^ b ^ T ' and ~dri U ^ 
where F % are the constraints ()2.2|) defining the m-phase solutions of (jl.lj) . The derivatives 
d/dk a , d/duj a , d/dn l here are the total derivatives including the dependence of the func- 
tions $ l (0, k, u}, n) on the corresponding parameters view the form of ()2.1()j) . However, all 
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these derivatives are identically equal to zero on the family A, so we get the first part of 
the Lemma. To finish the proof we have to note also that the values S(fc)(X, T), n^){X, T) 
do not appear in the k-th order of e of the normalization conditions ()2.11|) - ()2.12|) either. 

Lemma 2.1 is proved. 

We can formulate now the procedure in the following form: 

We try to solve the systems ()2.13|) recursively in all orders and find the functions 
X, T) satisfying the conditions (j2.11|) - ()2.12|) . At each k-th step of our procedure 
we get a system on the functions S(fc_i)(X, T), n^_i^(X, T) from the solvability conditions 
of ()2.13j) . The full set of Whitham solutions will then be parameterized by the set of all 
functions {S(fe)(X, T), n^(X, T)}, k > satisfying the conditions of solvability of systems 

dm. 

Now let us investigate the systems arising on the functions S(fc_i)(X, T), n^_i)(X, T) 
from the solvability conditions of (|2.13j) 

/2k i>2k _ <7 m /9 

<^ ]4 (*,X,T)jJ fc) (0,X,T) = 0, q =l,...,m + s (2.14) 

in the k-th order of e. We note first that the solvability conditions of (J2.13)) for k — 1 give 
a nonlinear system on the functions S^(X,T), n^(X,T). It's not difficult to prove the 
following Lemma: 

Lemma 2.2. 

The system arising for k = 1 coincides with the Whitham system for the functions 
S (0) (X,T) ; n (0) (X,T). 

Proof. 

Indeed, using Lemma 2.1 it's not difficult to see that the discrepancy f(i)(0,X, T) 
differs from i^(6,X, T) just by the terms containing derivatives 6 x and 6qt- However, 
according to Lemma 1.1 these terms do not affect the orthogonality conditions ()1.12j) 
which then coincide with orthogonality conditions for f(i)(#, X, T). 

Lemma 2.2 is proved. 

The Whitham system (orthogonality conditions) for the functions S(o)(X, T), n( )(X, T) 
can be written in the following general form: 

W® [S (0) , n (0) ] = A$ (S ( o)x, S (0) T, up)) Sforr + 

+ B a ] ( S (0)X, S(o)T,n (0 )) S(o)XT + C a ] ( S (0)X, S( )T,n (0 )) Sf Q)xx + 

+ G; 9) (S (0)x ,S (0)T ,n (0) ) n[ 0)T + (S (0) x, S (0) t, n (0) ) n\ 0)x = (2.15) 
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with some functions A a , , Ca , G\ , H\ , q = 1, . . . , m + s, a = 1, . . . , m, 
Z = 1, ...,s . 

Let us prove now the following Lemma about the systems on S(m, no,), k > 1. 
Lemma 2.3. 

T/ie orthogonality conditions of f(k+i) X, T) to the functions K [|| 0) no ] $™ e ^ e fol- 
lowing linear systems for the functions S^)(X, T), n^(X, T), k > 1: 



5SUX',T) S (^ X > T > dX dT + J J Sn i(X',T>) n (^ X > T > dX dT 



V (k) [S(o)> • • • , S(fc_i), ii(o), • • • , n (fe-i)] (X, T) 



In other words, the functions S^)(X,T), n^){X,T) satisfy the linearized Whitham 
system on the functions S(p){X,T), n^(X,T) with additional right-hand part depending 
on the functions S (0 ), . . ., S (fc _i), n (0) , . . ., n (fc _i). 

Proof. 

Let us look at the terms in fo. +1 )(0, X, T) which contain the functions S^){X,T), 
n {k) (X,T) : 

1) As we proved in Lemma 2.1 the functions f^JO, X,T) contain only the terms 

depending on S( ), n( ). Easy to see that the functions fJ k+1 J0,X,T) will then contain 
the terms 



/ / (1) V ; S? k ,(X',T') dX' dT' + / / ( ) V ; nUX',T) dX' dT' 
J J 6Sfa(X',T') WV ' ; J J 5n\ 0) (X',T>) ^ y ' 

according to the form of the first term in ([2. 10)1 . 

2) There are terms containing the functions Soa, and the function ^^(6, X, T). 
All such terms can be written in the form: 



6LL n U (X,T) , 5LL n U (X,T) 

b ^ [X ' 1 ] 6S? 0) (X',T') + ' } 6n\ 0) (X',T') 



& {1) (0,X,T) dX'dT' 



where L^ So , is the linear operator (jl.lOj) given by the linearization of the system (|2.2j) 
on the family of m-phase solutions. 
3) There will the terms of the form 
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r r S 2 FL n JO, X, T) 

- ... [Sq '" q1V 1 S? k JX',T) SUX",T") dX' dT' dX" dT" - ... 

J J SS^X'^SS^X^T") {k)y 

where FL n J0, X,T) is the left-hand part of the system (|2.2j) . 

However, the sum of all such terms is equal to zero since they all correspond to the 
expansion of the "shift of parameters" S^(X,T), n^(X,T) on the family A where we 
have F\0,X,T) = identically 

Let us look now at the terms in the orthogonality conditions 

/2n p2tt _ d m f) 

■■■J Q < Q0] ^,X,T) f{ k+1) (0,X,T) ^ = (2.16) 

containing the terms l)-2). 
We have identically 

/2n p2iv rl m R 
... k\$ v(0,X,T) LL nU (0,0',X,T) -— — = (2.17) 

on A (where L{ Sq >ao]j (0, 0', X, T) is the "core" of the operator L\ SoMj {X,T))- 
Easy to see then that the inner product of K^ q ' with the terms 2) is equal to 



2vr p2n 



JO 



6k,\$ JO.X.T) 8k\H v (0.X,T) 

b (k){X , 1 ) + n (k) (X , 1 ) 



X 



5S^(X',T) 
L\s ,^(X,T) ¥ {0,X,T) 



6n\ 0) (X>,T>) 



dX'dT' x 



d m 6 



It's not difficult to see now that the terms of orthogonality conditions containing the 



terms l)-2) can be written together in the form 
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5S? 0) (X',T') ^ [So ' no 



(9) 



f (1) [S 0) n ])(X,T) dX'dT' + 



Jo) 



(m) 



n \^ X ' X) 5n\ 0) (X>,T>) 



f (1) [S ,n ])(X,T) dX'dT 



(where <...-...> is the inner product of K^ q > and fm). 

All the other terms in the orthogonality conditions (|2.16|) are the smooth functionals 
of S( ), • • ., S(fc_i), rim), • ■ -5 n (fc-i)) so we get the statement of the Lemma. 

Lemma 2.3 is proved. 

Let us consider now the systems ()1.1|) satisfying the special non-degeneracy conditions. 
Namely, we will assume that the corresponding Whitham system (J2.15j) can be resolved 
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w.r.t. to the highest T-derivatives of the functions S^)(X,T), n^(X,T) and written in 
the "evolution" form: 



CO 

°{0)TT 



n 



(0)T 



M" o)0 (S ( o)x, S (0 )T,n ( o)) S (0)xx + ^ )/3 (S (0 )x,S(o)T,n (0 )) S {0)TX + 

Smvr,n (0) ) n p {Q)x , a = l,...,m (2.18) 



+ P(0)p (S(0) \ • ^fU)7'- 
r (0)/3 (S(O)X, O(0)T, 



'(0)XI 

S(o)T,n (0 )) n p 
After the introduction of the variables k 



Sm-iT, ri( )) S^ TX + 



+ R \o) P (S(o)x,S(o)T,n (0 )) .. (I1) V . 



l,...,s 



(2.19) 



system fl2~T^l) - (l2~TH|l in the form 



"(0) 



J (0)X ' ^(0) 



S(o)t we can wr ite the 



"(0)T 



a; 



(0)X 



uf 0)T = M^ )f3 (k ( o),w (0 ),n (0 )) fc (0)x + N^ )p (k (0) , w (0 ), n (0) ) w (0)jr + 



+ P {o)p ( k (o),^(o),n (0) ) 



rr 



(0)X 



n 



(0)T 



T (W ( k (o),^(o),n (0 )) fc (0 w + L( 0)/3 (k(o),W( ),n(o)) cjf n , v 4 



+ R \o)p ( k (o),^(o),n (0) ) 



n 



(0)X 



i.e. in the form p. 15)1 . 

We will consider now the hyperbolic systems ()2.2U|) . i.e. such that the matrix 



(2.20) 



^T( k (o),^(o),n 











M {0) 


N (0) 




T (0) 


L( ) 





has exactly N = 2m + s real eigen- values with N linearly independent real eigen- vectors. 
For hyperbolic systems ()2.20|) it's natural to consider the Cauchy problem with the smooth 
initial data k (0) (X,0), o; (0) (X,0), n (0) (X,0) (or S (0) (X,0), S (0 )t(A,0), n (0) (X,0)). The 
smooth solution of ()2.2()j) exists in general up to some finite time To until the breakdown 
occurs. So we can write the zero (global in X) approximation for the solution ()2.1()j) just 
in the time interval where we have a smooth solution of the Whitham system. Using 
Lemma 2.3 it's not difficult to prove then the following Lemma: 



Lemma 2.4 
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For non-degenerate hyperbolic Whitham system \2.20\) and the global solution S(o)(A, T), 
ri(o) (X, T) defined on the interval [0, To] the higher orders approximations in A2.1U\) are 
all defined for all X and T G [0, T ) and are parameterized by the initial values S^(X, 0), 
S (k)T (X,0), n (fe) (X,0). 3 

Proof. 

Indeed, as follows from Lemma 2.3 the functions S^)(X, T), n^)(X,T) are denned by 
the initial values S^(X, 0), S^)t(X, 0), n^)(X, 0) and can be found from the linear system 
using the characteristic directions of ()2.20|) (defined by S( )(X, T), n( )(X, T)) provided 
that all the smooth solutions S( )(A, T), . . . , S^-i)(X, T), n^(X,T) . . . , n( fc _!)(X, T) 
and ^ r (i)(0, X, T), . . . , \j>( fe _ 2 )(0, X, T) exist on the interval [0,T ). According to Lemma 
2.1 and Lemma 2.3 we can find then the functions ty^ k _^(0, X, T) which are the local 
expressions (in X and T) of S (0) (X,T), ... , S (fc _i)(X,T), n (0) (X,T) ... , n (A ._ 1) (X,T) 
and their derivatives. Using the induction we then finish the proof of the Lemma. 

Lemma 2.4 is proved. 

According to Lemma 2.4 we can formulate now the following statement: 

For the initial system (jl.lj) having the non-degenerate hyperbolic Whitham system 
()2.2()|) the corresponding Whitham solutions ()2.1()j) (or ()2.4|) - fl2.5j) ) are defined by the 
initial values S(A, 0, e), St(A, 0, e), n(X, 0, e) and exist in the time interval [0, T ) defined 
by the Whitham system (jQDp and the initial data S (0) (A, 0) = S(X, 0, 0), S (0) t(A, 0) = 
S T (X,0,0), and n (0) (A,0) = n(A,0,0). 

The deformation procedure. 

Let us note now that the series ()2.10j) (or ()2.4|) - (j2.5|0 give in fact the one-parametric 
formal solutions of (|l.lj) with a parameter e. Let us rewrite now the solutions ()2.10j) in 
the form which gives the concrete (formal) solution of (|l.lj) and is not connected with 
the additional one-parametric e-family including this given solution. We omit now the 
e-dependence of functions S(X, T, e), n(X,T,e) (or put formally e = 1) and say that the 
Whitham solution is defined now by functions S(X, T), n(X, T) determined by the initial 
values S(A, 0), St(A, 0), and n(A, 0). (Let us keep here the notations X and T for the 
spatial and time coordinates just to emphasize that we consider the "slow" functions 
S x (A,T),S T (A,T),n(A,T).) 

Thus, we define now the Whitham solution as the solution of the system 

F' l (ip, <p T , (fx, <Ptt, <Pxt, <Pxx, ■ ■ •) = , i = l,...,n (2.21) 
having the form 

3 Let us remind that we assume that all systems l|2.13|l are solvable if the corresponding orthogonality 
conditions are satisfied. 
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l (0,X,T) = ¥(S{X,T) + 0, S X (X,T), S T (X,T), n(X,T)) + 

+ ^$( fc) (S(X,T) + 0,X, T) (2.22) 

fe>i 

where the functions $| fc s : 

1) Are 27r-periodic with respect to each 6 a ; 

2) Have degree k (introduced below); 

3) Satisfy the normalization conditions (|2.11|) - ()2.12j) which will be written now in the 
form: 

/2tt p2ir n tj rn f) 
...jf ^% a (6, S x ,S T ,n) $l k) (e,X,T) — = (2.23) 

...jf 5>i,(0, S x ,S T ,n) 9\ k) (0,X,T) j-^ = (2.24) 
k > 1, (a = 1, . . . , m, / = 1, . . . , s). 

Let us introduce now the gradation used for the formal expansion ()2.22j) . Namely, for 
the systems (II having the non-degenerate hyperbolic Whitham system ()2.20|) we put the 
following gradation on the functions Sx{X,T), St{X,T), n(X,T) and their derivatives : 

1) The functions k a (X,T) = S X {X,T), uj a (X,T) = S$(X,T), and n l (X,T) have 
degree 0; 

2) Every differentiation with respect to X adds 1 to the degree of the function; 

3) The degree of the product of two functions having certain degrees is equal to the 
sum of their degrees. 

In other words, for the parameters U = (k, u>, n) we have the gradation rule of 
Dubrovin - Zhang type, i.e. 

All the functions /(U) have degree 0; 
The derivatives U kX have degree k; 

The degree of the product of functions having certain degrees is equal to the sum of 
their degrees. 

We put now the evolution conditions to the functions S(X,T), n(X,T) having the 
form : 

Stt = a (k) (k, w, n, k x , uj x , n x ,...) (2.25) 

k>l 
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u t = Yl *iU ( k > w ' n ' kx ' ux > nx ' • • • ) ( 2 - 26 ) 

k>l 

where rf^ are general polynomials in derivatives kx, wj, n^, kxx, ^xx, nil, • • • 
(with coefficients depending on (k, w,n)) having degree fc. 

We now substitute the expansion (J2.22j) in the system (|2.21|) and use the relations 
(|2.25j) - (j2.26|) to remove all the time derivatives of parameters (k, oj, n). After that we can 
divide the system (j2.21|) into the terms of certain degrees and try to find recursively all 
the terms $ (fc) (0,X,T) for k = 1, 2, . . . . 

It is easy to see again that for any X, T) we will have the linear system anal- 

ogous to UZ3}, (l2~T3D . i.e. 

L\s(x,T)Mx,T)]j $ (fc)(0>^> r ) = f(k){Qi X iT) (2.27) 
where f(fc)(0, X, T) is the discrepancy having degree k according to the definition above. 

We have to put again the orthogonality conditions 

■■■J Q ^{O^T) fafrX,!*) ^ = (2.28) 

on the functions ff k \(0, X,T) and then find the unique $(fy(0,X, T) satisfying the nor- 
malization conditions (J2~23|) - (j2~2l| . 

It's not difficult to prove the following Lemma: 

Lemma 2.5. 

1) For any system M.l)) having the non- degenerate hyperbolic Whitham system A2.2U\) 
the functions cr^, rf^ are uniquely determined by the orthogonality conditions \2.2ty in 
the order k. 

2) The functions a^, rjL-, give the Whitham system \2. 1^) - T2. 19\) for the functions 
S {0) (X,T), n (0) (X,T). 4 

Proof. 

Indeed, using Lemma 2.1 it is easy to see that the functions X,T) coincide 

with the functions f^(0, X,T) introduced in ()2.13|) after the replacement of functions 
S(o)(X, T), ri(o) (X, T) by S(X, T), n(X, T). Comparing then the orthogonality conditions 
(I2.28j) with (j2.14j) we get the second part of the Lemma. 

4 In fact, the functions S(X, T, e), n(X, T, e) introduced previously satisfy the full system (|2.25(1 - (|2.26|) . 
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To prove the first part we just note that all er^L, rf^ arise in the fe-th order of system 
(|2.21jl " in the same way" . We can conclude then that the orthogonality conditions (|2.28jl 
in the fc-th order always contain the functions cr^, rjl k s in one particular way which 
coincides with the appearance of <7?U, tL\ in the Whitham system arising for k — 1. From 
the definition of the non-degenerate Whitham system we now obtain the first part of the 
Lemma. 

Lemma 2.5 is proved. 

Definition 2.1. 

We call the system i2. 25}) - Hi. 2b)) or the equivalent system 

k'f = LC *x 

= a h ( k > w > n > k *> wx,nx,...) (2.29) 
n T = ^2 tfk) ( k ' w> n ' kx ' nx ' ' ' • ) 

fe>l 

i/ie deformation of the Whitham system \2.1§) - \2~Tfy) (or \2.2tty ). 

The functions k a (X,T), u a (X,T), n (X,T) are the "slow" functions of the variables 
x and t and the system (|2.29jl gives the analog of the "low-dispersion" expansion in our 
case. The asymptotic solutions ()2.22j) of the initial system (j2.21|) are parameterized by 
the asymptotic solutions k a (X, T), u a (X,T), n l (X,T) of the system ()2.29|) and arbitrary 
(constant) initial phases 9q. As follows from our considerations above the solutions (|2.22j) 
give all the "particular solutions" (|2.4j) - (J2.5|) . however, they do not contain the additional 
information about the one-parametric e-family given by (|2.4j) - (J2.5|) . 

Remark 1. 

Let us note that the full set of parameters of m-phase solutions of is given 

by k, oj, n, and 6q. However the functions 6 (X,T) do not present as the parameters of 
solutions (|2.22j) in this approach. This shows in fact that the introduction of the functions 
9q(X,T) do not give "new" formal solutions of and is responsible for the additional 
e-dependence of one-parametric families ()2.4j) - (j2.5|) . Here they are "absorbed" by the 
total phase S(X,T) connected with the "particular" formal solution of (Jl.lj) . 

Remark 2. 

In our consideration we fixed some functions k, oj, n) on the family A (for 

each (k, u>, n)) as having zero initial phases. However, the choice of the functions 
k, oj, n) is not unique. In particular, the natural change of the functions 
k, oj, n) on A can be written in the form 
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&(0, k, w, n) = $(6» + o (k, w, n), k, w, n) (2.30) 

where #q (k, u>, n) are arbitrary smooth functions. 

It's not difficult to see also that the system ()2.29|) depends on the choice of the functions 
k, u), n) in the high (k > 2) orders. 

Let us give here the definition given by B.A. Dubrovin and Y. Zhang ( I65j) and 
connected with the "equivalence" of different infinite (or finite) systems. 

Definition 2.2 (B.A. Dubrovin, Y. Zhang). 

Consider the system of the form 

U t = ^^)(U,U,,U tt ...) , u=l,...,N (2.31) 

k>\ 

for arbitrary parameters U v where all VK* are smooth functions polynomial in \Jx, Uxx, 
. . . and having degree k. We say that two systems \2.31\) are connected by the triviality 
transformation (or equivalent) if they are connected by the formal substitution 

u v = E^«( U > u *> u **> •••) 

fc>0 

where all UY k s are smooth functions polynomial in XJx, Uxi, . . . and having degree k. 

Let us say actually that the definition of B.A. Dubrovin and Y. Zhang is applied usually 
to the whole integrable hierarchies and plays the important role in the classification of 
integrable hierarchies having the form (j2.31j) . We will prove here the following Lemma: 

Lemma 2.6. 

The deformations of the Whitham system \2.29\) written for the initial functions 
<&(#, k, lo, n) and <&'(0, k, u), n) connected by the transformation A2.3U\) are equivalent 
in Dubrovin - Zhang sense. 

Proof. 

Let us prove first the following statement: 

For any transformation ()2.30|) there exists a unique change of functions S a (X), n l (X): 
s , a = S a _ 0* (k) Wj n) + J2 S? k) (k, w, n, k x , . . . ) 

k>l 

n 'i = n l + ^2 n U(K n, k x ,...) (2.32) 

fe>i 

such that: 

1) All S(fe), n(fe) are polynomial in derivatives of (k, oj, n) and have degree k; 
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2) For any solution (l2~2^|) O[ S ,n](0, X, T)) of (l2~2TJ) the functions 
$ /J (S'(X, T) + 0, S' x , S' T , n') satisfy the normalization conditions 

/2-7T /"27T n fl m fj 

■ J q ^<(S' + ^S^,S^n')0f Sin] (0,X,T) _ = (2.33) 

/27T /-27T n r\ m 

(2.34) 

For the proof of this statement let us note first that we can always express all the 
time derivatives of k, u>, and n using the system (|2.29|) in terms of X-derivatives of these 
functions. Using this procedure we try to find the transformation (|2.32|) recursively in 
all degrees by substitution of (|2.32|) in (|2.33|) - (j2.34|) . It's not difficult to check then that 
the functions Sfy, n l , k \ are defined in the k-th order of ()2.33|) - (j2.34|) from a linear system. 
The matrix of this linear system coincides with the Gram matrix of functions <lv , & t 
every X and T. Thus, for the full non-degenerate family of m-phase solutions of (jl.l)) this 
system has a unique solution at every degree k. The transformation 1)2.32)1 is evidently 
invertible in sense of the infinite series (polynomial w.r.t. derivatives of (k, to, n)) so we 
can also express the functions (S, n) in terms of (S', n'). 

We can try to use now the functions $'*(S' + 6, S' x , S' T , n') as the zero approximation 
in the (S', n')-expansion of the corresponding solution ()2.22j) . It is not difficult to see that 
the difference 

$"(S' + e,S^,S>') - ^ n] (9,X,T) 

can be represented as the infinite series polynomial w.r.t. derivatives of (k, oj, n) and 
starting with the terms of degree 1. After the expression of the functions (k, oj, n) in 
terms of (k', a/, n') in this difference we get finally the (S', n')-expansion of the solution 

The functions (S', n') satisfy now the new deformed Whitham system ()2.25|) - (j2.26|) 
corresponding to the choice of the functions $'(#, k', u/, n') as the functions of zero ap- 
proximation. 

Easy to see also that the transformation (|2.32j) remains polynomial in X-derivatives 
of (k, u), n) after the expression of time derivatives of (k, uj, n) using the system (|2.29|) . 

We obtain then that the transformation ()2.32|) gives a "triviality" connection between 
the systems f|2.29|) written for the initial functions <3?(0, k, u:, n) and &'(0, k, uj, n) . 

Lemma 2.6 is proved. 

At the end let us note again that the Lemma 2.6 is important in fact for the integrable 
hierarchies rather than for the one particular system (jl.l)) according to Dubrovin - Zhang 
approach to classification of integrable systems. 
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